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Abstract 

The  understanding  of  thin  film  growth  processes  and  their  control  requires  the  development  of 
surface-sensitive  real-time  optical  characterization  techniques  that  are  able  to  provide  insight  into 
the  surface  reaction  kinetics  during  an  organometallic  deposition  process.  These  insights  will  allow 
us  to  move  the  control  point  closer  to  the  point  where  the  growth  occurs,  which  in  a  chemical 
beam  epitaxy  process  is  a  surface  reaction  layer  (SRL),  built  up  of  physisorbed  and  chemisorbed 
precursor  fragments  between  the  ambient  and  film  interface.  This  contribution  presents  results 
on  parameter  estimations  of  rate  constants  and  optical  response  factors  in  a  reduced  order  surface 
kinetics  (ROSK)  model,  which  has  been  developed  to  describe  the  decomposition  and  growth  kinetics 
of  the  involved  organometallic  precursors  and  their  incorporation  in  the  film  deposition.  As  a 
real  time  characterization  technique,  we  applied  p-polarized  reflectance  spectroscopy  (PRS)  during 
low  temperature  growth  of  epitaxial  Gai_xInxP  heterostructures  on  Si(001)  substrates  by  pulsed 
chemical  beam  epitaxy  (PCBE).  The  high  surface  sensitivity  of  PRS  allows  us  to  follow  alterations 
in  composition  and  thickness  of  the  SRL  as  they  are  encountered  during  periodic  precursor  supply. 
The  linkage  of  the  PRS  response  to  the  ROSK  model  provides  the  base  for  the  parameter  estimation, 
giving  insights  into  the  organometallic  precursor  decomposition  and  growth  kinetics. 


*  Corresponding  author.  Tel.:  (404)463-9617;  Fax:  (404)651-1427;  E-mail:  ndietz@gsu.edu 


1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2000 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2000  to  00-00-2000 

4.  TITLE  AND  SUBTITLE 

Surface  Reaction  Kinetics  of  Gal-xInxP  Growth  During  Pulsed  Chemical 
Beam  Epitaxy 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

North  Carolina  State  University, Center  for  Research  in  Scientific 
Computation, Raleigh, NC, 27695-8205 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

The  original  document  contains  color  images. 

14.  ABSTRACT 

see  report 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

ABSTRACT 

18.  NUMBER 

OF  PAGES 

16 

19a.  NAME  OF 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


1.  Introduction 


The  engineering  of  advanced  optoelectric  devices,  microelectronic  circuits,  and  integrated  optical 
circuits  requires  precise  control  of  lateral  dimensions  and  thicknesses  of  device  features  and  perfection 
of  heterostructures.  Low  pressure  deposition  methods,  such  as  chemical  beam  epitaxy  (CBE)  and 
plasma  enhanced  chemical  vapor  deposition,  play  an  important  role  in  the  manufacturing  of  nanos¬ 
tructure  devices  and  advanced  ultra-large  scale  integration  (ULSI)  processing,  respectively.  Some 
areas  of  interest  include  improved  densities  of  integrated  electronic  devices  and  methods  of  improving 
the  control  of  epitaxial  deposition  to  realize  these  devices.  Key  targets  in  III-V  compound/silicon 
heterostructures  are  the  control  and  understanding  of  defect  formation  and  the  interactions  and 
propagations  of  defects  during  later  stages  of  compound  heteroepitaxy  growth.  To  understand  the 
properties  of  defects  in  the  epitaxial  film,  we  must  understand  the  kinetics  of  thin  film  growth. 

Understanding  and  controlling  thin  film  growth  has  been  a  slow  process  because  little  is  known 
about  chemical  reaction  properties  and  reaction  kinetics  parameters  during  the  decomposition  process 
of  metalorganic  precursors.  In  the  film  deposition  process,  the  film  surface  plays  a  major  role  in  the 
growth  process  as  the  addition  of  new  reactants  can  significantly  modify  the  growth  chemistries  [1,  2], 
This  characteristic  of  thin  film  growth  chemistry  during  chemical  deposition  combined  with  the 
demand  for  stringent  thickness  and  composition  control  in  advanced  optoelectronic  integrated  circuits 
has  led  to  the  development  of  surface  sensitive  real-time  optical  sensors  [3,  4]  that  are  able  to  move 
the  monitoring  and  control  point  closer  to  the  point  where  growth  occurs,  which  in  a  chemical  beam 
epitaxy  process  is  the  SRL.  The  SRL  is  built  up  of  physisorbed  and  chemisorbed  precursor  fragments 
between  the  ambient  and  film  interface. 

Applying  optical  probe  techniques  to  real-time  characterization  of  thin  film  growth  includes  the 
challenge  of  relating  surface  chemistry  processes  that  drive  the  growth  process  to  growth  and  film 
properties,  such  as  composition,  instantaneous  growth  rate,  or  structural  quality.  Most  characteri¬ 
zation  techniques  take  advantage  of  the  numerous  probes  available  to  provide  a  detailed  assessment 
of  the  ambient  by  accurately  measuring  ambient  process  parameters,  but  they  are  limited  in  their 
ability  to  deal  with  complex  surface  chemistry  processes.  For  example,  during  heteroepitaxial 
GaP/Gai_xInxP  growth  on  Si  under  pulsed  CBE  (PCBE)  conditions,  the  surface  is  exposed  to  a 
periodic  supply  of  metalorganic  precursors,  which  results  in  an  SRL  with  periodically  altered  thick¬ 
ness  and  composition.  The  measurement  technique  of  p-polarized  reflectance  spectroscopy  (PRS) 
[5,  6,  7],  which  has  demonstrated  high  sensitivity  towards  surface  reaction  processes  in  the  context 
of  real-time  monitoring  of  PCBE,  allows  us  to  undertake  the  desired  characterization  and  control  of 
thin  film  deposition. 

We  start  in  Section  2  with  a  brief  background  on  the  experimental  growth  and  monitoring 
conditions  and  show  results  obtained  by  PRS  during  real-time  characterization  of  heteroepitaxial 
growth  of  Gai_xInxP  on  Si  substrates.  In  Section  3,  we  introduce  the  model  used  to  simulate  the 
PRS  measurements.  We  also  describe  the  link  of  the  PR  response  to  the  simulation  parameters 
accessible  through  the  reduced  order  surface  kinetics  (ROSK)  model,  which  has  been  developed  to 
describe  the  decomposition  and  growth  kinetics  of  the  involved  organometallic  precursors  [4,  8].  The 
process  of  identifying  these  parameters  is  explained  in  Section  4.  Section  5  analyzes  the  results  of  the 
parameter  identification,  allowing  us  to  establish  and  validate  surface  reaction  kinetics  parameters, 
while  advancing  our  understanding  of  fundamental  chemistry  processes  in  thin  film  growth  using 
organometallic  precursors.  We  finish  with  our  concluding  remarks  in  Section  6. 
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2.  Experimental  Setup  and  Results 


For  monitoring  both  the  bulk  and  surface  properties  during  heteroepitaxial  Gai_xInxP  growth 
on  Si,  PRS  has  been  integrated  into  the  pulsed  chemical  beam  epitaxy  system  that  is  schematically 
shown  in  Figure  1.  In  PCBE,  the  surface  of  the  substrate  is  exposed  to  pulsed  ballistic  beams  of 


UHV  Chamber 


TMI 

Figure  1:  Setup  of  PCBE  system  for  III-V  compound  semiconductor  growth. 

TBP  [(C4H9)PH2],  TEG  [Ga(C2H5)3],  and  TMI  [In(CH3)3]  at  typically  350-450°  C  to  accomplish 
nucleation  and  overgrowth  of  the  silicon  by  an  epitaxial  GaP/Gai_xInxP  film.  For  PRS  and  laser 
light  scattering  (LLS),  we  employ  p-polarized  light  beams  with  wavelengths  A  =  632.8  nm  and  650  nm 
at  two  angles  of  incidence  (PR70:  4?i  =  71.5°  and  PR75:  4?i  =  75.8°).  Further  information  on  the 
experimental  conditions  is  provided  in  previous  publications  [4]- [9]. 

Each  growth  run  is  started  with  a  preconditioning  period,  where  the  substrate  temperature  is 
raised  to  the  growth  temperature  and  the  PR  signals  change  according  to  the  temperature  depen¬ 
dency  of  the  substrate  dielectric  function.  We  use  the  PR  signals  to  verify  independent  temperature 
measurements  and  calibrate  the  actual  surface  temperature.  A  constant  flow  of  Palladium  purified 
H2  (10  seem)  is  introduced  to  the  growth  chamber  during  the  preconditioning  and  growth  periods. 
The  background  pressure  is  10~9  torr  and  the  pressure  during  growth  is  in  the  range  of  10_3-10-4  torr. 

The  growth  sequence  starts  typically  after  1200  s  of  preconditioning  with  250  cycles  of  GaP 
growth  to  aid  in  lattice  matched  Gai_xInxP  growth  on  a  Si(001)  substrate.  After  the  deposition  of 
a  GaP  layer,  the  transition  to  Gai_xInxP  was  made  in  one  of  three  ways.  The  first  transition  was  an 
abrupt  change  to  the  desired  TMLTEG  flow  ratio.  The  other  two  transitions  are  graded  changes, 
done  linearly  either  manually  or  under  computer  control,  from  a  TMLTEG  flow  ratio  of  0  to  the 
desired  flow  ratio  (0-1).  This  flow  ratio  relates  to  the  composition,  x,  of  Gai_xInxP  as  discussed  in 
Section  3. 

Figure  2  shows  the  PR  and  LLS  signals  during  heteroepitaxial  growth  of  Gai_xInxP  on  Si(001). 
After  initiating  growth,  the  minima  and  maxima  observed  in  the  PR  signals  are  due  to  the  interference 
phenomena  as  the  film  grows.  Note  that  the  maxima  and  minima  of  the  two  signals  are  phase  shifted, 
due  to  the  fact  that  one  angle  of  incidence  (PR75)  is  above,  and  the  other  (PR70)  below  the  pseudo- 
Brewster  angle  of  the  growing  film.  Superimposed  on  the  interference  fringes  is  a  fine  structure, 
which  is  due  to  the  periodically  modulated  composition  and  thickness  of  the  SRL. 

For  the  surface  reaction  kinetics  analysis  and  the  validation  of  performed  simulations  using  the 
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Figure  2:  Monitoring  of  heteroepitaxial  Gai_xInxP  growth  under  PCBE  by  PRS  and  LLS.  The  inset 
enlargement  shows  the  PR  response  to  a  periodically  modulated  SRL  composition  and  thickness 
during  pulsed  precursor  supply. 


ROSK  model  presented  in  Section  3,  we  varied  one  experimental  parameter,  the  TMI  flow  rate. 
Each  growth  condition  was  set  to  the  specified  TMLTEG  flow  ratio  and  monitored  for  at  least  l\ 
interference  oscillations  in  order  to  get  stable  steady-state  growth  and  gather  sufficient  information 
to  analyze  and  simulate  the  growth  process. 

The  correlation  of  the  fine  structure  evolution  with  the  pulsing  sequence  of  the  precursor  supply 
is  depicted  in  Figure  3.  In  Figure  3,  the  PR  response  is  taken  during  steady-state  growth  on  a  rising 
flank  of  an  interference  oscillation  using  a  pulse  cycle  sequence  of  6  s,  a  TBP  pulse  from  0.0-0. 8  s 
and  3.0-3. 8  s,  a  TEG  pulse  from  1.5-1. 8  s,  a  TMI  pulse  from  4. 5-4. 8  s,  and  a  continuous  hydrogen 
flow  during  the  complete  sequence.  In  this  set  of  experiments,  the  flow  rates  and  pulse  durations 
of  TBP  (1.0  sccrn  for  800  ms)  and  TEG  (0.054  seem  for  300  ms)  were  kept  constant,  as  was  the 
pulse  duration  of  TMI  (300  ms).  The  flow  rate  for  TMI  was  varied  from  0.014  to  0.068  seem  to 
produce  the  desired  TMLTEG  flow  ratios  and  compositions,  x,  of  Gai_xInxP.  As  the  TMLTEG 
flow  ratio  increases,  the  magnitude  of  the  fine  structure  increases  according  to  a  larger  quantity  of 
TMI  constituents  in  the  SRL.  The  results  are  discussed  in  more  detail  in  Section  5,  which  includes 
comparisons  with  the  results  of  simulations. 
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Figure  3:  A  comparison  of  PR75  response  to  periodically  modulated  TBP,  TEG,  and  TMI  precursor 
pulses  for  different  TMLTEG  flow  ratios  (0.25-0.5-0.75-1.0-1.25)  at  the  rising  flank  of  an  interference 
oscillation. 


3.  Modeling  of  PRS  Results 


To  model  the  optical  responses,  the  growing  heteroepitaxial  film  is  approximated  by  a  four  layer 
model  consisting  of:  (1)  ambient  -  (2)  surface  reaction  layer  (SRL)  -  (3)  film  -  (4)  substrate.  We 
consider  here  Gai_xInxP  on  a  Si/GaP  virtual  substrate.  The  complex  reflectivity  coefficient  for 
p-polarized  incident  light,  given  a  four- layer  stack,  is  given  by  the  formula 


_  ri2 (1  +r-23r34e  2lf}a)  +  (r23  +  r34e  2^3)e  2l/?2 
r?  (1  +  r^^e-2*#5)  +  ri2(r-23  +  r34e^2^3)e^2lf32  ’ 

where  the  Fresnel  coefficients  rk(k+i)(k  =  1,2,3)  for  each  interface  are  given  by 


rk(k+ 1) 


£fc+i  V^k  ~  eisw24»i  -  efcy/efc+i  ~  eism24>i 
efc+iVefc  -  eism24>i  +  ekyjek+i  -  eisin2$i 


(1) 

(2) 


and  the  phase  angles  /3k  for  the  SRL  (k  =  2)  and  the  growing  film  epilayer  (k  =  3)  are  given  by 


Pk 


2  7T 

T 


dkVek~ 


(3) 
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In  equations  (l)-(3),  the  reflectivity  coefficient  rp  is  a  function  of  d2  and  d3  (the  thickness  of  the  SRL 
and  film  respectively),  ei,  €2,  63,  and  64  (the  complex  dielectric  functions  of  the  ambient,  SRL,  film, 
and  substrate  respectively),  $1  (the  angle  of  incidence),  and  A  (the  wavelength  of  the  impinging  laser 
light).  For  additional  information,  we  refer  the  reader  to  [9]. 

The  values  of  ei,  €3,  €4,  <f>i,  and  A  are  all  constant  in  time,  but  €2,  d2,  and  d3  vary  in  time  as 
the  film  grows  and  the  SRL  changes  in  composition  and  thickness.  To  understand  how  these  values 
change,  we  need  a  representative  model  of  the  chemical  kinetics  of  the  SRL  which  approximates 
the  pyrolysis  of  the  source  vapor  molecules  as  described  previously  [10,  11,  12].  For  TBP,  TEG, 
and  TMI  as  source  vapors  forming  Gai_xInxP,  we  employ  a  reduced  order  surface  kinetics  (ROSK) 
model  [4,  8].  The  ROSK  model  makes  the  following  simplifying  assumptions: 


•  TBP  pyrolysis  is  approximated  through  one  dominant  reaction  equation. 

•  TEG  decomposition  is  simplified  to  two  dominant  reaction  equations. 

•  The  formation  of  GaP  is  combined  into  one  dominant  reaction  equation. 

•  TMI  decomposition  is  simplified  to  two  dominant  reaction  equations. 

•  The  formation  of  InP  is  combined  into  one  dominant  reaction  equation. 

•  The  process  is  driven  by  a  periodic  source  vapor  supply  as  described  in  Section  2. 

Thus  the  kinetic  model  representing  the  SRL  reactions  can  be  represented  by  the  following 
system  of  ordinary  differential  equations: 


dt 

d 

dt 

d 

dt 

d 

dt 

d 

dt 

d 

dt 

d_ 

dt 


ni(t) 
n2{t) 
n3  (t) 
n4(t) 
n5(t) 
nett) 
n7(t) 


Si(t)  -  fcini(t)  -  kGapni(t)n3(t)  1Q_8mol  “  hnpni(t)nG(t)  1Q_8mol  (4) 

S2(t)  -  kdTEGn2(t)  -  k2n2(t)  (5) 

k2n2(t )  -  k3n3(t )  -  kGaPni(t)n3(t) j— y  (6) 

kGapni(t)n3(t)  1Q_8mol  (7) 

S 3  (t)  -  kdTMin5  (t)  -  k5n5(t)  (8) 

k3n3(t)  -  k6n6{t )  -  hupni(t)n6(t)  10-8mol  (9) 

fcinpm(t)n6(t)  1Q— 8mol  (10) 


The  units  of  each  term  of  these  equations  is  mol/s.  The  reaction  constants,  fcj,  have  units 
of  s-1,  while  the  molar  quantities,  rq,  have  units  of  mol.  The  state  variables  represent  the  number 
of  moles  of  the  various  components  of  the  SRL,  such  as  n\  for  active  surface  phosphorus  fragments. 
The  intermediate  gallium  precursor  fragments  in  the  TEG  decomposition,  possibly  diethylgallium 
(DEG),  are  represented  by  n2  and  the  active  gallium  fragments  are  represented  by  n3.  Similarly, 
the  intermediate  indium  precursor  fragments  and  the  active  indium  fragments  are  represented  by  n3 
and  uq  respectively.  In  equation  (4),  the  change  in  active  surface  phosphorus  fragments  is  written 
as  the  sum  of  a  source  term  Si,  a  desorption  loss  term  — k\n\ ,  and  reaction  terms  representing  the 
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formation  of  GaP  and  InP.  Equation  (5),  which  represents  the  defragmentation  process  of  TEG, 
contains  a  source  term  S2,  a  desorption  loss  term,  and  a  term  of  decomposition  from  the  intermediate 
gallium  precursor  fragments  to  active  surface  gallium  fragments.  Equation  (6),  which  represents 
the  change  active  surface  gallium  fragments,  has  a  term  of  creation  from  the  intermediate  gallium 
fragments,  a  desorption  loss  term,  and  a  reaction  term  for  the  formation  of  GaP.  Equation  (8), 
which  describes  the  defragmentation  of  TMI,  contains  a  source  term  S3,  a  desorption  loss  term,  and 
a  term  of  decomposition  from  the  intermediate  indium  precursor  fragments  to  active  surface  indium 
fragments.  Equation  (9),  which  represents  the  change  in  active  surface  indium  fragments,  has  a 
term  of  creation  from  the  intermediate  indium  fragments,  a  desorption  loss  term,  and  a  reaction 
term  for  the  formation  of  InP.  The  variables  714  and  717  represent  the  number  of  moles  of  created 
GaP  and  InP,  respectively,  integrated  in  the  deposited  Gai_xInxP  film  layer.  Equation  (7)  contains 
only  a  single  reaction  term  for  the  formation  of  GaP  from  active  surface  Ga  and  P,  and  equation  (10) 
contains  only  a  single  reaction  term  for  the  formation  of  InP  from  active  surface  In  and  P.  These 
two  equations  also  account  for  any  surface  activation  processes. 

The  source  terms  in  the  differential  equations  are  based  on  the  source  vapor  pulses.  These  are 
modeled  similarly  to  the  GaP  case  [13],  with  the  TBP  source  given  by  the  expression 


Si(t) 


Pi(£)7/3tbp 

Vtbp 


(11) 


where  Pi(t)  is  the  source  vapor  injection  rate,  Vtbp  is  the  molar  volume  of  TBP,  and  the  constant 
/3tbp  is  the  sticking  coefficient  of  TBP.  The  geometrical  parameter  7  represents  the  amount  of 
source  vapors  that  actually  hit  the  surface  of  the  wafer,  which  is  a  constant  value  for  a  given  reactor 
geometry.  Similarly,  the  second  and  third  source  terms  are  modeled  by 


S2(t) 

Ss(t) 


P2(*)7/3teg 

Vteg 

-P3(£)t/3tMI 

Vtmi 


(12) 

(13) 


with  corresponding  P2(t),  VteGj  and  /3teg  for  the  TEG  pulse,  and  P3(t),  Vtmi,  and  /3tmi  for  the 
TMI  pulse.  The  same  constant  7  appears  in  all  source  terms.  For  each  source  term,  we  use  a 
constant  flow  rate  between  start  and  stop  times,  and  zero  flow  elsewhere.  Due  to  the  time  needed 
for  the  source  vapor  gates  to  open  and  close  and  for  the  vapors  to  travel  to  the  surface,  there  is  a 
small  time  difference  between  the  start  (and  stop)  of  the  vapor  pulse  and  the  start  (and  stop)  of 
the  source  vapor  arrival  at  the  surface.  We  account  for  this  time  difference  with  a  parameter  delay. 
With  a  source  vapor  pulse  starting  at  ton  and  stopping  at  t0g,  the  source  vapors  actually  reach  the 
surface  starting  at  t0 n  +  delay  and  ending  at  tQff  +  delay.  We  estimated  a  delay  of  0.27  s  using  the 
parameter  estimation  process  as  described  in  Section  4. 

The  system  of  differential  equations  (4)-(10),  together  with  the  source  terms  ( 1 1)-(  13) ,  and 
appropriate  initial  conditions,  can  be  solved  numerically  for  the  number  of  moles  ri\,  n2,  n 3,  714,  715, 
Tig,  and  717.  From  these  solutions,  the  film  and  SRL  thicknesses  are  obtained  from  the  formulas  [4] 

d3{t)  =  -j[VGa,prH{t)  +  Einpn7(t)]5  (14) 

d2(t)  =  ^p[Vi7ii(t)  +  V2n2{t)  +  V3n3(t )  +  V5n5(t )  +  V6n6(t)},  (15) 

where  A  is  the  surface  area  of  the  wafer,  14  are  the  molar  volumes  of  the  components  7i&,  and  EcaP 
and  Vjnp  are  the  molar  volumes  of  GaP  and  InP,  respectively.  Finally,  agjjx  is  an  effective  SRL 
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thickness  parameter  representing  the  percentage  of  the  SRL  that  contributes  to  the  optical  response. 
The  effective  dielectric  function  of  the  SRL  is  expressed  as  the  sum  over  all  molar  fractions  x,-L 
contributing  to  the  SRL  and  is  given  by 


e2{t)  =  1  +  ^2  xi{t)Fi  and  xi(t) 
4,7 


nj(t ) 

E k  nk{ty 


(16) 


which  is  derived  from  the  Sellmeier  equation  [4,  14]. 

With  the  values  of  the  temporal  dependent  parameters  e2l  d2 ,  and  d3  found  by  equations 
(14)-(16)  and  with  the  constant  parameters  ei,  €3,  €4,  4>i,  and  A,  the  reflectivity  coefficient  rp  can 
be  calculated  using  equations  (l)-(3).  From  rp.  we  then  find  the  value  that  is  actually  measured  in 
the  experiments  by  computing  the  reflectance  Rp  =  \rp\2. 

Finally,  the  composition,  x,  for  the  compound  semiconductor  Gai_xInxP  is  expressed  as  the 
averaged  molar  ratio  of  molar  concentration  over  a  cycle  sequence  [4]: 


rt—6  sec  d  ,r.  JJ. 

Jt= 0  sec 


X  = 


ft=  0  sec  (ln4  +  in7)dt 


The  instantaneous  film  growth  rate,  gji,  is  given  by 

9fl 


1 

A 


T  _  d  d 

LGaP^-^4  +  Finp— n7 


(17) 


(18) 


In  the  equations  described  above,  4>i,  A,  eu  Vi,  V2,  V3,  V5,  V6,  VGap,  RinP,  Rtbp,  Lteg,  Rtmi,  A 
/3tbp;  /3teG;  /3tmIj  «srl;  and  all  start /stop  times  and  flow  rates  contributing  to  Pi,  P2,  and  P3  are 
known  quantities.  The  values  of  the  dielectric  functions  63  and  €4,  the  rate  constants  k\,  k2.  k3, 
fc.5,  fcg,  fegTEG)  ^dTME  ^GaP,  and  fcinp,  the  optical  responses  Pi,  F2.  P3,  P5,  and  Fq,  the  geometrical 
parameter  7,  and  delay  have  to  be  determined.  In  the  next  section,  an  inverse  least  squares  problem 
will  be  formulated  to  estimate  these  parameters.  The  comparison  of  the  model  using  the  estimated 
parameters  with  the  experimental  data  will  be  done  in  Section  5. 


4.  Three  Layer  Problem  and  Setup  of  Parameter  Identification 


We  formulate  here  an  inverse  least  squares  problem  for  identifying  the  unknown  parameters  in 
the  model  in  Section  3.  We  want  the  set  of  parameters  that  results  in  the  modeled  reflectance  most 
closely  matching  the  experimental  data.  For  this  we  look  for  the  vector  of  parameters 

Q  =  {Fi,F2,F3,F5,F6,ki,  fcdTEG  A2A3,  ^GaP ,  &dTMi  A5,h,  hnP ,  7 ,  delay)  (19) 

that  minimizes  the  cost  function 


m 


I  E(*  exp  (ti)  -  Rca\c{ti,q)Y 


(20) 


In  this  cost  function,  Rexp(ti)  is  the  experimental  reflectance  data  set  at  the  measurement  times  ti, 
and  RcaicYi^O)  =  \rp |2,  where  rp  is  the  complex  reflectivity  coefficient  given  by  (1),  is  the  simulated 
reflectance  calculated  at  the  same  times  using  the  parameter  set  q. 
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Since  larger  numbers  of  parameters  make  the  minimization  process  increasingly  difficult,  we  do 
not  include  €3  and  e 4  in  q.  We  are  able  to  remove  these  two  parameters  from  our  large  parameter 
estimation  problem  by  formulating  a  seperate  but  simpler  estimation  problem.  The  simpler  problem 
uses  a  three-layer  stack  as  a  model  of  the  growing  film  by  removing  the  SRL  from  consideration  and 
leaving  just  the  ambient,  film,  and  substrate  layers.  This  allows  us  to  first  study  the  properties 
of  the  film  itself  before  later  looking  at  the  behavior  in  the  SRL.  The  formula  for  calculating  the 
reflectance  for  a  three-layer  stack  is  given  by 

_  n 3  +  r34e~2^3 

3,P  1  +  ?'l3?'34e_2*^3  ’ 

where  r  13  and  ?’34  are  Fresnel  coefficients  for  the  reflection  from  interfaces  1-3  and  3-4  (now  that 
layer  2,  the  SRL,  is  removed),  and  the  phase  angle  /33  is  for  the  film  layer.  These  values  are 
calculated  by  formulas  analogous  to  (2)  and  (3). 

To  compare  results  from  this  formula  with  experimental  results,  we  use  a  method  described 
previously  in  our  investigation  with  GaP  formation  [13].  At  this  point  we  focus  on  the  large- 
amplitude  interference  oscillations,  which  have  a  periodicity  of  several  hundreds  of  seconds  and 
are  formed  by  the  film  growth,  and  remove  the  small-amplitude  fine  structure  oscillations  that 
are  modulated  with  the  precursor  injection.  After  extracting  the  experimental  three-layer  stack 
reflectance  i?3,exP)  we  can  identify  the  parameters  e3  and  €4,  as  well  as  the  average  growth  rate,  gr, 
which  is  used  to  find  the  film  thickness  at  times  U  ,  by  comparing  the  calculated  reflectance  Alcaic  = 
]?’3)Pj2  from  equation  (21)  to  i?3,exp-  This  is  done  through  an  inverse  least  squares  formulation  by 
finding  r  =  (e3,  €4,  gr)  that  minimizes  the  cost  function 

J(r)  =  (-R3,exp(ii)  -  ^calc^i,^))2-  (22) 

The  results  of  the  three-layer  curve  fitting  provide  us  with  the  average  growth  rate,  gr ,  and  the 
film  dielectric,  e3,  for  each  TMLTEG  flow  ratio.  The  substrate  dielectric,  €4,  is  also  found,  but 
is  not  flow  ratio  specific  since  it  takes  into  account  the  initial  GaP  layer  and  the  abrupt  or  graded 
change  in  precursor  pulse  flow  to  the  desired  TMLTEG  flow  ratio  of  Gai_xInxP  film  growth.  Specific 
parameter  estimates  for  different  flow  ratios  are  presented  in  Table  1. 

025  05  075  TO 

0.31  0.41  0.61  0.80 

10.56-1. 30i  10.94-1.48i  10.53-1. 08i  10.15-1. 50i 

Table  1:  Average  values  of  flow  ratio  specific  three-layer  model  parameters 

Once  we  have  determined  the  values  of  €3,  €4,  and  gr,  these  parameters  are  used  to  solve  the 
four-layer  stack  parameter  identification  problem,  which  finds  the  unknown  parameters  F\,  F%,  F$,  F5, 
Fq,  k\ ,  A’d'nx;,  ^2,  £'3,  ^GaP,  fcdTMl,  h,  kG,kinp,-f,  and  delay. 

5.  Gai_xInxP  Growth  Kinetics  Analysis 

We  compare  measurements  taken  with  various  TMLTEG  flow  ratios  while  all  other  conditions 
are  fixed  so  that  we  may  determine  the  dielectric  function  (€2)  corresponding  to  a  given  compostion,  x, 


1.25 

093 

1 1 .28-0.651 


TMLTEG  flow  ratio 
gr(A/a) 

£3 
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of  Gai_xInxP.  In  addition,  varying  the  TMLTEG  flow  ratio  reveals  several  important  characteristics 
in  the  fine  structure.  We  will  explain  these  features  and  show  how  the  mathematical  model  of  the 
growth  process,  using  the  reduced  order  surface  kinetics  model  described  in  Section  3,  captures  these 
phenomena. 

Looking  at  the  fine  structure  in  the  PR  75  data,  the  most  noticeable  change  that  occurs  with 
the  TMLTEG  flow  ratio  variation  is  the  increasing  of  the  fine  structure  amplitude  as  TMLTEG  flow 
ratio  increases,  as  illustrated  in  Figure  4.  The  regions  near  the  TBP  pulse  exhibit  an  upward  slope  in 


Figure  4:  Experimental  and  simulated  PR75  responses  for  various  TMLTEG  flow  rates  under  steady- 
state  growth  conditions. 

the  fine  structure,  while  the  regions  near  the  TEG  pulse  exhibit  a  downward  slope,  as  was  described 
in  work  with  GaP  [13].  The  upward  slope  in  the  optical  response  is  due  to  the  source  TBP  exposure, 
the  defragmentation,  and  the  generation  of  active  phosphorus  on  the  surface,  while  the  downward 
slope  is  due  to  the  source  TEG  exposure,  the  defragmentation,  and  the  active  gallium  generation  on 
the  surface.  The  regions  in  the  optical  response  fine  structure  related  to  a  TMI  exposure  exhibit  a 
downward  slope,  much  like  the  TEG  exposure,  except  that  the  deflection  is  on  a  much  larger  scale 
for  the  larger  TMI  flow  rates. 

In  Figure  4  we  can  also  examine  the  changes  in  the  simulated  and  experimental  data  as  we  varied 
the  TMLTEG  flow  ratio.  Here  we  see  that  the  amplitude  of  the  fine  structure  increases,  especially 
after  the  TMI  pulse,  as  we  increase  the  flow  ratio.  In  addition,  the  steeper  slope  in  the  interference 
fringes  is  related  to  the  increase  in  the  overall  growth  rate  due  to  the  greater  TMI  flow.  Both  the 
experimental  and  simulated  data  sets  exhibit  these  characteristics  and  the  results  agree  with  each 
other  very  well. 
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One  larger-scale  feature  observed  in  the  reflectance  data  is  change  in  the  average  film  growth 
rate  for  the  various  TMLTEG  flow  ratios.  The  frequency  of  the  large-scale  interference  fringes 
is  also  related  to  the  TMLTEG  flow  ratio  and  the  average  film  growth  rate.  Equations  (2),  (3), 
and  (21)  illustrate  that  a  larger  growth  rate  corresponds  to  a  more  quickly  increasing  dj,  and  /?3, 
which  increases  the  frequency  of  the  interference  fringes. 

The  other  large-scale  characteristic  feature  of  the  data  set  is  the  position  of  (or  the  lack  of) 
turning  points  in  the  fine  structure,  which  can  be  depicted  in  the  first  derivative  of  the  reflectance 
data,  as  shown  in  Figure  5.  The  turning  point  positions  are  characterized  by  the  PR  measurements 


Time  (s) 


Figure  5:  The  first  derivatives  of  the  simulated  and  experimental  reflectance  data,  shown  with  the 
experimental  data  set. 

through  the  places  where  the  amplitude  of  the  reflectance  derivative  is  minimized.  There  will  either 
be  pairs  of  turning  points  for  each  interference  oscillation,  as  the  fine  structure  switches  sign,  or 
a  single  minimum  in  the  case  that  there  is  no  sign  change.  Since  the  turning  point  positions  are 
functions  of  the  SRL  dielectric  properties,  the  good  match  between  the  derivatives  in  Figure  5  implies 
that  the  SRL  dielectric  properties  have  been  well  represented  by  the  model. 

The  steps  in  the  generation  of  the  set  of  simulated  data  for  flow  ratio  0.75  presented  in  Figure  4 
are  shown  in  detail  in  Figure  6.  The  five  SRL  components  are  the  result  of  the  source  pulses  and  the 
ROSK  simulations.  Once  we  have  the  SRL  components,  we  find  the  SRL  thickness  and  dielectric 
function.  These  values  then  are  used  to  calculate  the  reflectance.  Figure  6  shows  how  the  exposure 
of  the  growth  surface  to  TBP  changes  the  SRL  and  with  it  the  optical  response,  causing  an  upward 
slope  in  the  fine  structure.  Similarly  the  TEG  and  TMI  exposures  cause  changes  in  the  SRL  which 
result  in  downward  slopes  in  the  optical  response.  The  good  fit  of  the  simulated  fine  structure  to  the 
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Time  (s) 

Figure  6:  The  steps  in  the  generation  of  a  set  of  simulated  data  are  each  represented  by  their 
contributing  numerical  solutions.  In  order  from  top  to  bottom:  (i)  the  flux  of  precursor  pulses,  (ii) 
the  molar  quantities,  rq(t),  (iii)  the  SRL  thickness,  e?2(i),  (iv)  the  SRL  dielectric  function,  62  (t),  (v) 
the  simulated  and  experimental  PR  responses. 


experimental  data  shown  in  Figure  6  is  also  obtained  along  the  entire  interference  oscillations.  The 
agreement  of  the  model  with  the  experimental  data  in  the  several  ways  discussed  above  supports  the 
efficiency  of  the  ROSK  model  for  describing  the  growth  process  and  the  observed  optical  response. 

The  parameters  used  for  the  minimization  process  in  the  determination  of  the  unknown  param¬ 
eters  in  the  model  are  presented  in  Table  2.  The  geometrical  parameter  7  =  0.024  was  estimated  in 
the  minimization  process  for  a  two  inch  diameter  circular  substrate  wafer. 

The  numerical  solutions  were  calculated  using  programs  written  in  MATLAB  code.  The  stiff 
system  of  differential  equations  in  the  ROSK  model  was  solved  numerically  by  the  built-in  function 
“ode23s”,  an  adaptive  mesh  and  low-order  Runge-Kutta  method  [15].  The  least  squares  optimization 
problem  was  solved  using  a  Nelder-Mead  [16,  17,  18]  algorithm. 
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molar  volumes  (cm6 /mol) 

sticking  coefficients 

effective  SRL  thickness 

Vtbp  =  128.57 

V'i  =  17.02 

Ptbp  =  0.15 

otsRL  =  0.75 

Vteg  =  148.36 

V2  =  13.00 

Pteg  =  1-0 

Vtmi  =  101.99 

C3  =  11.80 

Ptmi  =  1-0 

VcaP  =  12.21 

V5  =  13.00 

VInP  =  15.22 

V6  =  15.71 

Table  2:  Constants  used  in  the  calculations  of  the  model 


We  estimated  the  following  parameters  in  the  four  layer  stack  model  by  averaging  the  re¬ 
sults  of  independent  best  fits  of  the  experimental  data  sets  corresponding  to  different  TMI:TEG 
flow  ratios:  the  rate  constants  k\  =  3.57  s-1,  k^TEG  =  0.60  s-1,  k2  =  1-87  s-1,  k 3  =  0.05  s-1, 
k GaP  =  1-86  s-\  kdTMi  =  0.52  s-1,  k§  =  0.62  s-1,  =  0.52  s_1,  and  kjnp  =  1.58  s-1;  and  the 
optical  responses  F\  =  12.77  —  0.05i,  F2  =  13.09  —  0.92i,  F$  =  7.96  —  4.66i,  F$  =  11.43  —  l.lli,  and 
Fq  =  4.52  —  9.93i.  The  averages  of  the  SRL  dielectric  function  for  these  parameters  are  given  in 
Table  3  for  each  TMI:TEG  flow  ratio.  Figure  7  shows  the  film  composition  from  the  model  plotted 
against  experimental  data  as  a  function  of  the  TMLTEG  flow  ratios. 


TMLTEG  flow  ratio 

0.25 

0.5 

0.75 

1.0 

1.25 

C2 

12.32-1. 51i 

12.29-l.55i 

12.39-1.38i 

12.21-l.63i 

11.99-1.90i 

Table  3:  Average  values  of  the  SRL  dielectric  function 


Figure  7:  Film  compositions  for  the  simulations  using  various  TMLTEG  flow  ratios,  as  well  as  some 
experimental  data  points. 
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6.  Conclusions 


We  extended  our  previous  effort  on  utilizing  a  reduced  order  surface  kinetics  model  using  gen¬ 
eralized  reaction  rate  parameters  for  the  growth  of  GaP  [13]  to  describe  the  decomposition  kinetics 
of  the  organometallic  precursors  TBP,  TEG,  and  TMI  used  during  the  heteroepitaxial  growth  of 
Gai_xInxP  on  Si/GaP.  A  set  of  coupled  differential  equations  was  used  to  describe  the  surface  re¬ 
action  kinetics  of  our  three  precursor  pulses  and  provided  information  about  the  dynamics  of  molar 
concentrations  of  precursor  fragments  stored  in  the  surface  reaction  layer  and  their  incorporation  into 
the  growing  film.  The  model  was  validated  and  the  unknown  parameters  involved  were  identified 
by  fitting  simulated  PRS  measurements  from  the  model  with  several  experimental  data  sets.  The 
results  demonstrate  that  a  mathematical  model  can  be  effectively  used  to  describe  the  film  growth 
process.  Current  efforts  are  directed  towards  molecular  specific  optical  process  monitoring  to  vali¬ 
date  the  predicted  surface  reaction  layer  constituents  and  their  concentrations  as  computed  by  the 
ROSK  model. 
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